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We study the kinetic freeze-out conditions of hadrons in 
Au+Au collisions at 11.6 AGeV/c using different parametri- 
zations of an expanding thermal fireball. We take into ac- 
count the available double differential momentum spectra of 
a variety of particle species, covering a large fraction of the 
total momentum space. The overall fit to the data is very 
good and indicates a relatively low kinetic freeze-out temper- 
ature of about 90 MeV with an average transverse expansion 
velocity at midrapidity of about 0.5 c. 

PACS numbers: 25.75.-q, 25.75.Ld, 24.10.Pa 

When studying the nuclear phase diagram in extreme 
regions of temperature and density with heavy-ion colli- 
sions one is faced with the difficulty that, due to the li- 
mited size and strong dynamics of the fireballs created in 
such collisions, full (global) thermodynamic equilibrium 
and the thermodynamic (infinite volume) limit can never 
be reached. The best one can hope for is a state of lo- 
cal thermal equilibrium in a hydrodynamically expanding 
environment, with thermodynamic homogeneity volumes 
of several 100 to 1000 fm'^ [|l[. If sufficient local equili- 
bration is achieved one can try to extrapolate from the 
experimental findings to the thermodynamic limit. The 
most exciting prospect of such an endeavour is the exper- 
imental confirmation of a deconfinement phase transition 
in hot and dense nuclear matter 

In order to reach this goal an extensive experimental 
program has been launched studying heavy-ion collisions 
at various beam energies. The onset of a new phase of 
nuclear matter would probably be most clearly seen in 
a study of different observables as functions of beam en- 
ergy. For this purpose it is necessary to carefully analyze 
the data at all available beam energies. Here we present 
a study of the experimental particle spectra from Au-|- Au 
collision at a beam momentum of 11.6 A GeV/c measured 
by several groups at the Brookhaven AGS. The goal 
is to find a simple, but realistic parametrization of the 
freeze-out stage in these collisions. The extracted freeze- 
out parameters, especially the temperature and the mean 
transverse fiow velocity, will be compared with other col- 
lision systems and with collisions at different energies. 
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This should help to understand the gross features of the 
collision dynamics at ultrarelativistic energies. 

For the description of particle production we start from 
the formalism of Cooper and Frye Q which describes the 
single-particle spectrum as an integral over a freeze-out 
hypersurface, thus summing the contributions from all 
space-time points at which the particles decouple from 
the fireball: 

Here g is the degeneracy factor and f{x,p) the momen- 
tum distribution at space-time point x. In thermal mod- 
els one uses for f{x,p) a thermal equilibrium distribution 
and determines Ef by a freeze-out criterium for thermal 
decoupling Q . For the low temperatures discussed below 
the Boltzmann approximation is sufficient, but we allow 
for a space-time dependence of the temperature T, the 
chemical potential /i, and the flow velocity u^^: 

/(:.,p)=exp(^ j . (2) 

Our choice of the geometry of the freeze-out surface 
Sf and of the space-time dependence of the parameters 
reflects a compromise between relative simplicity of the 
model and importance of its dynamical ingredients. Fo- 
cussing on central collisions, we assume azimuthal sym- 
metry of the spatial geometry and momentum distribu- 
tions. We further assume boost-invariant collective dy- 
namics along the longitudinal {z) direction This was 
shown before ^ to give a better representation of the 
measured rapidity spectra than a static "Hagedorn" flre- 
ball. Also, a static fireball yields a much stronger ra- 
pidity dependence of the transverse slope parameters, 
Tsiopo — r/cosh(y — ?/mid), than observed The as- 

sumption of longitudinal Bjorken flow suggests to use lon- 
gitudinal proper time r = ^/t"^ — and space-time rapid- 
ity 77 = Artanh (z/i) as suitable variables in the t-z-plane; 
the transverse radial coordinate is denoted by r . Midra- 
pidity is defined to be at 77 = 0. Longitudinal boost- 
invariance of the flow fleld is implemented by the choice 
u^(x) — '^(X^v xXx^er^Vzix)) with Vz{T,r±,ri) = tanhr;. 
In transverse direction we take a linear flow rapidity pro- 
flle, tanhwj_ = ^(^7)^, where Rq is the transverse radius 
at midrapidity and p{ri) will be specified below. 

The geometry of the freeze-out hypersurface Sf is fixed 
as follows: In the time direction we take a surface of 
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constant proper time, r = tq. In ry-direction the freeze- 
out volume extends only to a maximum space-time ra- 
pidity 77inax; this breaks longitudinal boost-invariance, as 
required by the finite available total energy. In the trans- 
verse direction the boundary is given by R{ri) which may 
depend on rj. We consider two choices: 



R{7^) = i?o • e (77; 
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The first choice describes a cylindrical fireball in the 
ry-r _L-space. For 200 A GeV S-l-S collisions at the SPS 
such a picture was shown to work well |Q, but for the 
lower AGS beam energy this can no longer be expected. 
An elliptic or cigar shape as given in Eq. should be 
more appropriate |lO| . 

Having specified the freeze-out geometry and the dis- 
tribution function (||) at freeze-out we obtain the follow- 
ing thermal single-particle spectrum: 
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Before comparing with data we must still add the con- 
tributions from resonance decays. This is done following 
Refs. [|ll|,|l2|, again using thermal distributions (||) for 
the resonances. This procedure implies the assumption 
of full chemical equilibrium, including strange particles. 
As we will see this assumption is not very well justified, 
but sufficient for estimating resonance feeddown contri- 
butions. We only include resonances with masses up to 
the A(1232) resonance. Tests showed that at the low 
freeze-out temperatures discussed below the inclusion of 
higher resonances has very little infiuence on the results. 

At this point only T{x) and ^i{x) remain to be spe- 
cified. The chemical potential of a hadron with baryon 
number B and strangeness S is fixed by the baryonic 
and strange chemical potentials as iJi{x) = B ■ ^^{x) -\- S ■ 
Hsix) (we neglect small isospin asymmetry effects @). 
We have investigated two scenarios. In the first we set 
for simplicity T, fiB, and fig constant along the freeze- 
out hypersurface. In the second we select a Gaussian 
77-dependence of the freeze-out temperature, 



T(?7) = Toe"(^) 



(6) 



and fix iJ-Biv) by requiring freeze-out at constant baryon 
density. This is reasonable because in the range of T and 
fiB discussed below meson-meson reactions occur much 



less frequently than meson-baryon or baryon-baryon re- 
actions, making the baryon density the relevant con- 
trol parameter for freeze-out. The strangeness chemi- 
cal potential is obtained as a function of T(r]) and 
fJ-B{v) from the condition of local strangeness neutral- 
ity ijl^. The thermal fit parameters in the model are 
thus ^b(j]—^) si-nd To, the baryon chemical potential and 
temperature at midrapidity, and the width A77 of the 
temperature profile (g). 




FIG. 1. Proton and pion rriT-spectra from E866 The 
proton spectra are given in rapidity bins of width 0.1; the 
top spectra correspond to 0.5<j/iab<0.6, the bottom ones to 
l-6<j/iab<1.7. (Midrapidity is at i/iab = 1.6.) Successive spec- 
tra are scaled down by a factor 10. The pion spectra are 
given for 1.4<yiab<1.6; the vr" spectrum is divided by 10. 
Solid lines: best fit with Model 3. 

Boost-invariant longitudinal flow combined with a 
cylindrical fireball geometry leads to nearly rapidity- 
independent transverse mass spectra near midrapid- 
ity. (For the present data midrapidity corresponds to 
Z/iab=l-6.) This contradicts the measurements which 
show a clear rapidity dependence of the average trans- 
verse momentum (px) and of the transverse slope pa- 
rameters also near midrapidity, especially for the heavier 
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protons and A's |^-^ (see also Figs. and the discus- 
sion below). To describe the data one must therefore 
either reduce the transverse flow or the temperature as 
one moves away from midrapidity. Having investigated 
a wide range of possibilities, we here report on the re- 
sults obtained by combining in various ways the follow- 
ing relatively simple options: Cylindrical vs. ellip- 
tic (^ transverse geometry; constant vs. 77-dependent 
(^ freeze-out temperature; and a constant {p{r])=po) vs. 
an 77-dependent (p(77)=po-\/l— (?7^/?7max)) transverse flow 
gradient. The corresponding model parameters are listed 
in Table |. 



symbol description 

V ~ toRq fireball volume 

77max maximal longitudinal flow rapidity 

po maximal transverse flow rapidity 

A'b(7? ~ 0) baryon chemical potential at midrapidity 

To temperature at midrapidity 

Al] width of the longitudinal temperature profile 



TABLE L Model parameters. Note that tq and _Ro can be 
scaled out from the integrals in (|^) , leaving only the combina- 
tion ToRq as an independent fit parameter. It is proportional 
to the Lorentz invariant proper fireball volume ^ = Jj.. u^^dof^ 
which is quoted in Table 

Model 1 combines cylindrical geometry with constant 
temperature and transverse flow; as already mentioned 
it does not produce a very good fit. The parameters cor- 
responding to the smallest x^/DOF are listed in the first 
column of Tabic ||. The quoted value of xVDOF=6.64 
was obtained after allowing for different overall normal- 
izations for the various particle species than predicted 
by chemical equilibrium fl^ . The corresponding normal- 
ization factors Mi are also given; their deviation from 1 
indicates the departure from chemical equilibrium at the 
point of thermal freeze-out. Effects of up to 25% are seen 
and should be expected in view of the finding in 
that chemical freeze-out at the AGS occurs already at 
Tehom ^ 130 MeV. 

Model 1 fails in two characteristic ways: (1) The ra- 
pidity spectra for heavy particles are too flat near midra- 
pidity and too steep near y = ±77max- (2) The slopes 
of the TTiT-spectra are essentially constant in the region 
|y| ^ ?7max while the data show a significant steepening of 
the spectra already closer to midrapidity. Both failures 
are due to the boost-invariance of Model 1 in the region 

hi < ?7max- 

Due to the well-known weak sensitivity of the rapidity 
spectra to changes in the temperature |Q| it is hard to 
improve on them by simply allowing the freeze-out tem- 
perature to depend on 77; with the parametrization (|^) 
we found it to be impossible. A smooth reduction of the 
effective tranverse flow as described above works much 
better, in particular if the simultaneous improvement on 



the rapidity dependence of the transverse slopes is taken 



into account. 


Nevertheless, 


the proton and A 


rapidity 


distributions are still a bit too flat in the middle and too 










parameter 


Model 1 


Model 2 


Model 3 


l//(10''fm^) 


1.74 


1.74 


1.82 


^max 


1.05 


1.66 


1.71 


po 


0.726 


0.793 


0.792 


fisiv = 0)/MeV 


548 


549 


536 


To /MeV 


90.5 


91.6 


93.3 


Ari 


00 


00 


3.58 


/is (77 = 0)/MeV 


61.1 


62.7 


59.5 


{vi_){v = 0) 


0.449 


0.484 


0.483 




0.924 


0.920 


0.920 


K- 


1.101 


1.090 


1.078 




1.218 


1.239 


1.268 




0.887 


0.915 


0.914 




0.972 


0.978 


0.972 


A/'a 


1.115 


1.078 


1.206 


xVdof 


6.64 


3.16 


2.40 



TABLE II. Fit results for the different models (see text for 
their description). 



Surprisingly, a very good fit of both the rapidity and 
all the TOT-spectra can be obtained by simply reducing 
the transverse size of the flreball according to (^) , with- 
out changing the freeze-out temperature or transverse 
flow gradient {Model 2). The reduced transverse radius 
in backward and forward rapidity regions not only leads 
to a faster and smoother decrease of the proton and A 
rapidity distributions, but the implied reduced average 
transverse flow away from midrapidity allows for an ex- 
cellent fit to the steepening TJix-spectra. Combining Eq. 

with an 77-dependence of the transverse flow gradi- 
ent overdoes it: the rapidity spectra become too sharply 
peaked near midrapidity. Some further fine-tuning of 
Model 2 is possible by combining the elliptic geometry (^) 
with an additional 77-dependence of the freeze-out tem- 
perature [Model 3). However, the temperature gradients 
resulting from the fit are small, i.e. the width A7/ 3> ?7max 
is large. In spite of the considerable expense in comput- 
ing time resulting from the introduction of the additional 
parameter A77 its effect on x^/DOF is only minor. 

In Figures 1-3 we compare the results of Model 3 with 
the data; the corresponding curves from Model 2 would 
have been hardly distinguishable. Remaining minor flaws 
of Model 3 are a slightly too broad rapidity spectrum 
and somewhat too flat mx-spectra in the very backward 
region for protons and a slight lack of curvature in the 
TTiT-spectra for pions. Additional pion 77iT-spectra cov- 
ering the backward rapidity region 0.6 < j/iab < 1-4 have 
recently become available |T7|; although we did not in- 
clude them in our fit they agree with the predictions of 
our Model 3. Indeed, the lack of concavity of our calcu- 
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lated pion mx-spectra becomes less of a problem at back- 
ward rapidities. 




"T "'0 '"T 

FIG. 2. Kaon mx-spectra from E866 M, in rapidity bins 
of width 0.2. Successive spectra are scaled down by a factor 
10. The top spectra correspond to 0.6<yiab<0.8, the bottom 
ones to 1.8<yiab<2.0. Solid lines: best fit with Model 3. 



70 
60 
50 
:o40 
?30 
20 
10 

50 

:o4o 

?30 
20 
10 




f \ ^ 

k \ 
1 \ 




/ \ 


f \ 



10 



-1 



60 
50 
40 
30 
20 
10 




12-2-1012 

y-Vmid y-ymid 

FIG. 3. Rapidity spectra of protons, A hyperons and pions 
in Au-|-Au collisions at the AGS. The A's were measured b 
E891 1^. For protons and pions we use data from E866 
near midrapidity and from E877 for forward rapidities. 
The measured data points were reflected about midrapidity 
using the symmetry of the collision system. Solid lines: best 
fit with Model 3. 



parametrization . Our analysis thus strongly suggests 
an average kinetic freeze-out temperature of about 90 
MeV and an average radial flow velocity at midrapidity 
close to 0.5 c. Similarly low kinetic freeze-out tempera- 
tures were found in jl^] for the smaller Si-|-Au system, 
in a model which has many similarities with ours. How- 
ever, the analysis of |0 also included information from 
two-particle Bose-Einstein correlations which were not 
available for our system. The average transverse flow ve- 
locity at midrapidity in Si-|-Au is a bit smaller: with a 
linear transverse velocity profile Chapman and Nix pT| ] 
found (vj_){ri = 0) = I X (0.683 ± 0.048) c « 0.45 c. 

In a successful fit of the rapidity spectra for pions, 
kaons, protons and A's with a freeze-out temperature of 
130 MeV was reported. The rapidity spectra are, how- 
ever, not very sensitive to the freeze-out temperature 0; 
we found that the large body of mx-spectra are badly 
fit with T = 130 MeV and clearly point to a consider- 
ably lower kinetic freeze-out temperature. On the other 
hand, the thermal freeze-out parameters found here cor- 
respond to a hadron resonance gas with an entropy per 
baryon S/A ~ 12—13. This is the same range as found 
in the chemical freeze-out analysis of Si-|-Au collisions at 
the AGS |]l5| and confirmed later in Refs. |16 ^,|l^ also 
for Au-|-Au collisions. Our findings are therefore con- 
sistent with chemical freeze-out at Tchem ~ 130 MeV, 
followed by isentropic (hydrodynamic) expansion and ki- 
netic freeze-out at Tkin ~ 90 MeV. 

Naively one would expect that increasing beam energy 
leads to higher energy deposition in the reaction zone and 
thus to higher temperature and/or stronger radial flow 
at freeze-out. A tendency for stronger flow at larger ^/s 
is indeed seen at energies below and including the AGS 
|^,|l|. At the higher SPS energy an analysis of nega- 
tive particle spectra and two-particle correlations from 
the NA49 Pb-|-Pb experiment points to kinetic freeze- 
out at T=100-120 MeV with (u_l)=0.45-0.5 c 
the lower freeze-out temperature and larger flow being 
preferred This supports our expectation, but con- 

tradicts an earlier analysis 1^,^ which did not differ- 
entiate between chemical and kinetic freeze-out and thus 
gave larger freeze-out temperatures and smaller flow ve- 
locities. 

We thank J. P. Wessels and C.A. Ogilvie for kindly 
sending us their data files, and J. Rafelski for a very con- 
structive remark. We gratefully acknowledge the support 
of this work by BMBF, DFG and GSI. 



It is interesting to note that the thermal parame- 
ters, including the average radial flow velocity are 
quite stable against the discussed variations in the model 
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